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Abstract 

An (n, d) permutation array (PA) is a subset of S n with the property 
that the distance (under some metric) between any two permutations in 
the array is at least d. They became popular recently for communication 
over power lines. Motivated by an application to flash memories, in this 
paper the metric used is the Chebyshev metric. A number of different 
constructions are given as well as bounds on the size of such PA. 



1 Introduction 

Let S n denote the set of all permutations of length n. A permutation array of 
length n is a subset of S n . Recently, Jiang et. al [HE] showed an interesting new 
application of permutation arrays for flash memories, where they used different 
distance metrics to investigate efficient rewriting schemes. Under the multi- 
level flash memory model, we find the Chebyshev metric very appropriate for 
studying the recharging and error correcting issues. We note it by d max . For 
7r, a G S n , d m ax( 7r i ') = niaxj |7r j; — <jj|. We consider a noisy channel where 
pulse amplitude modulation (PAM) is used with different amplitude levels for 
each permutation symbol. The noise in the channel is an independent Gaussian 
distribution with zero mean for each position. The received sequence is the 
original permutation distorted by Gaussian noise, and its ranking can be seen 
as a permutation, which can be different from the original one. 

To study the correlations between ranks, several metrics on permutations 
were introduced, such as the Hamming distance, the minimum number of trans- 
positions taking one permutation to another, etc. [3J, [!]• For instance, Stoll 
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and Kurz [5] investigated a detection scheme of permutation arrays using Spear- 
man's rank correlation. Chadwick and Kurz [6] studied the permutation arrays 
based on Kendall's tau. 

Under the model of additive white Gaussian noise (AWGN) [7] , there is only 
a small probability for any amplitude level to deviate significantly from the 
original one. This inspired us to use the Chebyshev distance. Observe that 
two permutations with a large Hamming distance can actually have a small 
Chebyshev distance and vice versa. They appear to complement each other in 
some sense. 

In this paper, we give a number of constructions of PAs. For some we 
give efficient decoding algorithms. We also consider encoding from vectors into 
permutations. 

2 Notations 

We use [n] to denote the set {1, . . . , n}. S n denotes the set of all permutations 
of [n]. For any set X, X n denotes the set of all n-tuples with elements from X. 

Let l denote the identity permutation in S n . The Chebyshev distance be- 
tween two permutations n, a € S n is 

dmax(7T, O") = max{|7Tj - (Jj | | 1 < j < n} . 

An (n, d) permutation array (PA) is a subset of S n with the property that 
the Chebyshev distance between any two distinct permutations in the array is 
at least d. We sometimes refer to the elements of a PA as code words. 

The maximal size of an (n, d) PA is denoted by P(n, d). Let V(n, d) denote 
the number of permutations in S n within Chebyshev distance d of the identity 
permutation. Since d max (''7 c) = d max (ir, ira), the number of permutations in S n 
within Chebyshev distance d of any permutation tt 6 S n will also be V(n,d). 
Bounds on P(n, a) and V(n, d) will be considered in Sec. |4j 

3 Constructions 

In this section we give a number of constructions of PAs, one explicit and some 
recursive. 

3.1 An explicit construction 

Let n and d be given. Define 

C = {(7Ti, . . . , 7r n ) E S n \^i = i (mod d) for all i £ [rt]}. 
If n — ad + b, where < b < d, then C is an (n, d) PA and 

|C| = {{a+l)\) b (a\) d - b . 
In particular, we get the following bound. 
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Theorem 1. Ifn = ad + b, where < b < d, then 

P(n,d) > ((a + l)l) b (a\) d - b . 
Example 1. For d = 2, we get 

P(2a,2) > (a!) 2 . 

We note that if 2d > n, then a = 1 and b = n - d and so \C\ = 2 n ~ d . If 
2d = n, then a = 2, b = 0, and we have C = 2 d = 2"~ d as well. However, if 
2d < n, then |C| > 2™~ d . Especially, when d is small relative to n, C is much 
larger than 2"~ d . For example, for n = 30, d = 2, |C|/2"- d « 6.37 x 10 15 . 

This construction has a very simple decoding algorithm. For d > 2t + 1, we 
can correct error up to size £ in any coordinate. For coordinate i, the codeword 
has value ixi = i (mod d). Suppose that this coordinate is changed into a = 
TTi + u, where \u\ < t. Then iri is the integer congruent to i which is closest to 
a. Therefore, decoding of position i is done by first computing 

a = i — a (mod d), 

where — (d — l)/2 < a < (d — l)/2. Then a = —u, and so we decode into 
a + a = TTi . 

3.2 First recursive construction 

Let C be an (n, d) PA of size M, and let r > 2 be an integer. We define an 
(rn, rd) PA, C r , of size M r as follows: for each multi-set of r code words from 

C, 

(7r^,...,7r«),j=0,l,...,r-l, 

let 

Pi = ( r7r i j) -.?>•••) r7T n ) - j), i = °: 1> • • • . r - !> 

and include (po|pi| • ■ ■ |Pr-i) as a codeword in C r . It is clear that under this 
construction the distance between any two distinct pj, py is at least rd. It 
is also easy to check that (po|pi| • • • |/V-i) G S rn . Hence \C r \ = P(n,d) r . In 
particular, we get the following bound. 

Theorem 2. If n > d and r > 2, then 

P(rn,rd) > P(n,d) r . 

3.3 Second recursive construction 

For a permutation n = (m, 7T2, . . . , 7r n ) € 5„ and an integer m, 1 < m < n + 1 
define 

^ m (7r) = (m,7ri,7T2, • • .,n' n ) G S^+i 
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by 

= %i if ■n l < m, 

-k[ = 7Tj + 1 if 7T. t > m. 

Let C be an (n, d) PA, and let 

1 < si < s 2 < • ■ ■ < s t < n + 1 

be integers. Define 

C[si,S2, ...,s t ] = {ip Sj {n)\l<j<t 7 ne C}. 
Theorem 3. // C is an (n, d) PA of size M and 

sj + d < Sj + i for 1 < j < t — 1, 
then C[si, S2, . . . , s t ] is an (n + 1, d) PA of size tM . 

Theorem 4. // C is an (n, d) PA of size M and n < 2d, then C[d] is an 
(n + I, d + 1) PA of size M. 

Proof. If j > j', then 

d max ((p Sj (it), ip Sj , (o-)) > Sj - Sj> > d. 

Next, consider j' — j. If n, a G C, it =/= er, then w.l.o.g, there exist an i such 
that TTi > o~i + d. Hence 

, , , , / TTi - (Ti + 1 > d if 7Tj > Sj > Oi, 

d m ^ Sj (tt), ^ W) > ( ^ _ ^ > d otherwise. 
This proves Theorem [3] To complete the proof of Theorem 0] we note that 

TTi > o~i + d > d + 1 > d, 

and 

o~i < TTi — d < n — d < d. 

Hence 7r.; > c? > <Ti and so 

"max 

{<f Sj (ir),tp Sj (a)) > d+1. 



The constructions imply bounds on P(n,d). First, choosing t = \n/d\ + 1, 
s t =n + l and Sj = (j — 1) [n/d\ + 1 for 1 < j < t — 1, we get the following 
bound. 



Theorem 5. If n> d> 1, then 

n > ( 

-d 



P(n + l,d) > ( ^ +l)P(n,d). 
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Example 2. In Example[l}we showed that the explicit construction implied that 
P{2a,2) > (a!) 2 . Combining Theorem\^and search, we can improve this bound. 
We have found that P(7, 2) > 582, see the table at the end of the next section. 
From repeated use of Theorem^ we get 

Q7 

P(2a,2)>( a ( a -l)...5) 2 .4F(7 ; 2)>-(a!) 2 . 

Theorem [4] implies the following bound. 

Theorem 6. If d < n < 2d, then 

P(n+l,d+l) > P(n,d). 

Theorem |S] shows in particular that for a fixed r, 

P(d + 1 + r, d + 1) > P(d + r, d) for d > r. (1) 

We will show that P(d + r, d) is bounded. We show the following theorem. 

Theorem 7. For fixed r, there exist constants c r and d r such that P(d-\-r, d) = 
c r for d > d r . Moreover, 

c r <2 2r {2r)\ (2) 

and 

d r < 1 + (2r - l)c. r - r. (3) 

Remark. The main point of Theorem [7] is the existence of c r and d r . The 
actual bounds given are probably quite weak in general. For example, Theorem 
[7] gives the bounds c\ < 8 and d\ < 8. In Theorem [5] below, we will show that 
ci = 3 and d\ — 2. Theorem [Jj gives C2 < 384 and < 1151, whereas numerical 
computation indicate that C2 = 9 and di = 5. 

We split the proof of Theorem [JJ into three lemma. 

Lemma 1. Ifd>r, then P(d + r,d) < 2 2r (2r)!. 

Proof. Suppose that there exists an (d + r,d) PA C of size M > 2 2r (2r)L 
We call the integers 

1, 2, . . . , r and d + 1, d + 2, . . . ,d + r 

potent, the first r smaller potent, the last r larger potent. Two potent integers 
are called equipotent if both are smaller potent or both are larger potent. If the 
distance between two permutations (tti,it2, . . . , 7r„), (p\, p%, . . . , p n ) is at least 
d, then there exists some position i such that, w.l.o.g, — pi > d, Then m 
is a larger potent element and pi is smaller potent. Each permutation in Sd+ r 
contains 2r potent elements and we call the set of positions of these the potency 
support xC 71 ") °f the permutation, that is, the potency support of 7r is 

x (tt) = {i | 1 < 7T j; < r} U {i \ d + 1 < tt,; < d + r}. 
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The potency support of C is the union of the potency support of the permuta- 
tions in C, that is 

x(C) ={i | 1 < 7T; < r for some ir G C} 

U {i | d + 1 < TTi < d + r for some tt G C}. 

Let 7r G C. For each p G C, p ^ 7r, we have <i(7r, p) > d Hence there exists 
some i G xl 71 ") such that pi is potent. Therefore, the set 

{(p,i) p G C and i G x(tt)} 

contains at least 2r + (M — 1) > M elements. Hence there is an i G x( 7r ) such 
that 

\{p £ C \ pi is potent}| > M/(2r) > 2 2r (2r - 1)!. 

Since 

{p G C | pi is potent} = {p G C | pi is smaller potent} 

U{p G C | p-i is larger potent}, 

there exists a subset Ci C C such that 

|d | >2 2r - 1 (2r-l)! 

and the elements in position are equipotent. 

We can now repeat the procedure. Let 7r G C\. There must exist an ii G 
xM \ {*i} such that 

{p G Ci | p l2 is potent}| > |C 1 |/(2r - 1) > 2 2r ~ 1 (2r - 2)!. 

Hence we get subset Ci C Ci such that 

|C 2 | > 2 2r ~ 2 (2r-2)! 

and the elements in position ii are equipotent (and the elements in position i\ 
are equipotent). 

Repeated use of the same argument will produce for each j, 1 < j < 2r a 
set Cj such that 

|Q| >2 2 ^(2r-j)! 

and for j positions i\, Z2, . . . ij, the elements in those positions are all equipotent. 
In particular, | C 2l - 1 > 1, all permutations in Cir have the same potency support 
{ii, %2, ■ ■ ■ , «2r}: and for each of these positions, all the elements in that position 
are equipotent. This is a contradiction since the distance between two such 
permutations must be less than d. Hence the assumption that a PA of size 
larger than 2 2r (2r)! exists leads to a contradiction. 

Lemma [1] combined with (fTJ) proves the existence of c r and d r and gives the 
bound ©. 
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Lemma 2. If C is a (d + r, d) PA of size M where 

d> r and d + r > \x(C)\, 

then there exists a (d—l+r, d—1) PA of size M . In particular, if M = P(d+r, d), 
then 

P(d-l + r,d-l) = P{d + r,d). 

Proof. Replace all elements in range r + 1, r + 2, . . . , d in the permutations 
of C by a star * which will denote "unspecified". The permutations in C is 
transformed into vectors containing the potent elements and d — r stars. Note 
that if we replace the unspecified elements in each vector by the integers r + 
1, r + 2, . . . , d in some order, we get a permutation, and the distance between 
two such permutations will be at least d since we have not changed the potent 
elements. 

Since the length d + r of C is larger than |x(C)|, there exists a position 
where all the vectors contains a star. Remove this position from each vector 
and reduce all the larger potent elements by one. This given a set of M vectors 
of length d—l + r and such that the distance between any two is at least d—1. 
Replacing the d — 1 — r stars in each vector by r + 1, r + 1, . . . , d — 1 in some 
order, we get an (d — 1 + r, d — 1) PA of size M. 

If M = P(d + r,d), then we get 

P(d - 1 + r, d - 1) > P(d + r,d). 
Since P(d - 1 + r, d - 1) < P(d + r, d) by ^ , the lemma follows. 

Lemma 3. // C is a (d + r, d) PA of size M and d> r, then 

\ X (C)\ <M(2r-l) + l. 

Proof. Each permutation has potency support of size 2r. The potency 
support of any two permutations in C must overlap since their distance is at 
least d. Hence each permutation after the first will contribute at most 2r — 1 
new elements to the total potency support. Therefore, 

|X(C)| <2r+(M-l)(2r-l). 



Remark. By a more involved analysis, we can improve this bound somewhat. 
For example, we see that two new permutations can contribute at most 4r — 3 
to the total support. 

We can now complete the proof of Theorem [7] Let C be a (d + r, r) code 
of size c r . By Lemma EJ |x(C)| < c r (2r - 1) + 1. If d > 1 + c r (2r - 1) - r, 
then d + r > |x(C)|. Hence, by Lemma [2j P(d -l+r,d-l) = P(d + r,d). 
Therefore, d r < 1 + c r (2r — 1) — r, that is, ([3]) is satisfied. This completes the 
proof of Theorem [7] 
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Theorem 8. We have P(d + 1, d) = 3 for d > 2. 

Proof. We use the same notation as in the proof of Lemma [2l Let C be an 
(d + l,d) PA. The only potent elements are 1 and n. W.lo.g. we may assume 
the first permutation in C is (1, n, .) where * denotes some unspecified 

integer in the range 2,3, ... ,d. W.l.o.g, a second permutation has one of three 
forms: 

(n,l, *,*,...), (n,*,l, *,...), (*, *,...)• 

We see that if the second permutation is of the first form, there cannot be more 
permutations. If the second permutation is of the form (n, *, 1, *, . . .), then there 
is only one possible form for a third permutation, namely (1, *, n, *, . . .). Hence 
we see that P(d + 1, d) < 3 and that P(d + 1, d) = 3 for d > 2. 

To determine P(d + r, d) along the same lines for r > 2 seems to be difficult 
because of the many cases that have to be considered. Even to determine 
P(d + 2, d) will involve a large number of cases. For example for the second 
permutation there are 138 essentially different possibilities for the four positions 
in the potency support of the first permutation. For each of these there are many 
possible third permutations, etc. 

3.4 Encoding/decoding of some PA constructed by the 
second recursive construction 

Suppose we start with the PA 

C d = {(1,2,3,..., d)}. 
For v = d, d + 1, . . . , n — 1 let 

C v+1 = C v [l,v + l]. 

Then C n is an (n, d) PA of size 2 n ~ d . For some applications, we may want 
to map a set of binary vectors to a permutation array. One algorithm for 
mapping a binary vector (x\, x-i, . . . , x n -d) into C„ would be to use the recursive 
construction of C n by mapping (x\,X2, ■ ■ ■ ,Xi) into a permutation 7r in Cd+i- 
Recursively, we can then map (x\, x-i, ■ ■ ■ , Xi, 0) to /i(7r) and (x\,xi, . . . , Xi, 1) 
to / d+i+ i(7r). 

However, there is an alternative algorithm which requires less work. Retrac- 
ing the steps of the construction, we see that given some initial part of length 
less than n — d of a permutation in C n , there are exactly two possibilities for the 
next element, one "larger" and one "smaller". More precisely, induction shows 
that if the initial part of length i — 1 contains exactly t "smaller" elements, then 
element number i is either t + 1 (the "smaller") or n — i + t + 1 (the "larger"). 
This is the basis for a simple mapping from Z^~ d to C n . We give this algorithm 
in Figure Q] 

We see that the difference between the larger and the smaller element in 
position i < n — d is n — i. Hence we can recover from any error of size less than 
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Input: (x u ...,x n - d ) £Z 2 " 

Output: (7T1, . . . , 7T„) G C„ 

for i <— n — d + 1 to n do Xj <— 0; 
i <— 0; //* t is the number of zeros seen so far.*// 
for z * — 1 to n do 
if iEj = 

then {-Ki <- t + 1; i <— t + 1;} 

else {7Tj <— n — i + t + 1;} 



Figure 1: Algorithm mapping Z^ 1 to C, 



Input: (tti, . . . ,n n ) £ [n] n 
Output: (xi, . . .,x n - d ) 
t <— 0; //* t is number of zeros determined. */ / 
for i <— 1 to n — d do 

if 7T, < (n - i)/2 + t + 1 
then {xi <— 0; t <— f + 1;} 
else {xi <— 1;} 



Figure 2: Decoding algorithm recovering the binary preimage from a corrupted 
permutation in C n . 
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(n—i)/2 by choosing the closest of the two possible values, and the corresponding 
binary value. We give the decoding algorithm in Figure [2j 

Without going into all details, we see that we can get a similar mapping 
from q-aiy vectors. Now we start with the PA 

c (? _ 1)d = {(1,2,3,..., te-rjd)}. 

For (q - l)d < v < n - 1 let Sj = (j - 1) \y/{q - 1)J + 1 for 1 < j < q - 1 and 
Sq = v + 1. Let 

C*i/+1 = C V \S\, S2, ■ ■ ■ , Sg\. 

Then C n is an (n, d) PA of size q n -(i~ 1 ) d . Encoding and decoding correcting 
errors of size at most {d— 1)/2, based on the recursion, is again relatively simple. 

4 Further bounds on P(n, d) 
4.1 General bounds 

Since d max (ir,<7) < n — 1 for any two distinct permutations in S n , we have 
P(n, n) = 1. Therefore, we only consider d < n. 

Since the spheres of radius d in S n all have size V(n, d), we can get a Gilbert 
type lower bound on P(n,d). 

Theorem 9. For n > d > 2 we have 

P(n,d) > 



V(n,d-1) 

Proof. It is clear that the following greedy algorithm produces a permutation 
array with cardinality at least n\/V(n,d — 1). 

1. Start with any permutation in S n . 

2. Choose a permutation whose distance is at least d to all previous chosen 
permutations. 

3. Repeat step 2 as long as such a permutation exists. 

Let C be the permutation array produced by the above greedy algorithm. Once 
the algorithm stops, S n will be covered by the |C| spheres of radius d— 1 centered 
at the code words in C. Thus nl < \P\ ■ V(n, d — 1) which implies our claim. 

Similarly, since the spheres V(n, [(d— 1)/2J) are disjoint, we get the following 
Hamming type upper bound. 

Theorem 10. Ifn>d>l, then 

P{n, d) < 



V(n, L(d-1)/2J) 
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If n < 2d and d is even, we can combine the bound in Theorem [TU] with 
Theorem [6] to get the following bound which is stronger than the ordinary 
Hamming bound, at least in the cases we have tested. 

Theorem 11. If d is even and 2d > n > d > 2, then 

Pin, a) < — ; tt^t- 

1 ' ' ~ V(n + l,d/2) 

Example 3. For n = 11 and d = 6, Theorem 1 1 01 gives 

p(n,6) < ,. = — r - = :«<;r 

whereas Theorem \11\ gives 

P(ll,6)< t— 1 — = . = .s:,(). 



11! 




11! 


_y(n,2)_ 




_11854_ 


12! 




12! 


Lv(12,3)J 




563172 



Remark. We can of course use Theorem [6] repeatedly r times and then 
Theorem HUl to get 



P{n,d) < 



(n + r)\ 



V(n + r, l(d + r- 1)/2J) 



for all r > 0. However, it appears we get the best bounds for r = 1 when d is 
even and r = when d is odd. 

In general, no simple expression of V(n, d) is known. A survey of known 
results as well as a number of new results on V(n, d) were given by Kl0ve [5]. 
Here we briefly give some main results. 

As observed by Lehmer [9], V(n,d) can be expressed as a permanent. The 
permanent of an n x n matrix A is defined by 



perA = ^ ' ' ' a n^„ 



In particular, if A is a (0, l)-matrix, then 

per A = \{n G S n : = 1 for all 

Let yl(" !<i ) be the n x n matrix with a-"' d ' = 1 if \i — j\ < d and a)"'^ = 
otherwise. 

Lemma 4. V(ra,d) = perA^™^). 
Proof. 

V(n,d) = \{tt £ S n ■ 7r) < d} 

= {7T G S'n : |i — 7Tj| < (i for all i}\ 
= |{tt G S n : a^f = 1 for all *}| 
= perA^"'*. 
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For fixed d, V(n, d) satisfies a linear recurrence in n. A proof is given in fTU] 
(Proposition 4.7.8 on page 246). For 1 < d < 3 these recurrences were deter- 
mined explicitly by Lehmer [5], and for 4 < d < 6 by Kfove 8 . In particular, 
this implies that 

lim V(n 1 d) 1 / n =fx d , 

n — >oo 

where fid is the largest root of the minimal polynomial corresponding to the 
linear recurrence of V(n,d). Lehmer [5] determined fid approximately for d — 
1, 2, 3 and Kfove [8] for d < 8. 

For an n x n (0, l)-matrix it is known (see Theorem 11.5 in |llj ) that 

n 

perA < H^;!) 1 ^, 

i=l 

where is the number of ones in row i. 

For we clearly have < 2d + 1 for all i. Hence 

V(n, d) < [(2d + l)l]«/( 2d + 1 ) for all n (4) 

and 

^<[(2d+l)!] 1/(2d+1) . 
In Table [T] we give [id and this upper bound. 



Table 1: fid and its upper bound. 



d 


[id 


[(2d + l)!]Vta*fi) 


Hd/(2d + l) 


1 


1.61803 


1.81712 


0.53934 


2 


2.33355 


2.60517 


0.46671 


3 


3.06177 


3.38002 


0.43739 


4 


3.79352 


4.14717 


0.42150 


5 


4.52677 


4.90924 


0.41152 


6 


5.26082 


5.66769 


0.40468 


7 


5.99534 


6.42342 


0.39969 


8 


6.73016 


7.17704 


0.39589 



We note that for large d, fid/ (2d + 1) w 1/e. 
Combining Theorem [9] and ((4|) we get 

Corollary 1. For n > d > 1, we have 
P(n,d) > 



[(2d- l)!]n/(M-l) • 
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4.2 Table of bounds on P(n,d) 

We have used the following greedy algorithm to find an (n, d) PA C: Let the 
identity permutation in S n be the first permutation in C. For any set of per- 
mutations chosen, choose as the next permutation in C the lexicographically 
next permutation in S n with distance at least d to the chosen permutations in 
C if such a permutation exists. The size of the resulting PA is of course a lower 
bound on P(n, d). 

The lower bounds in Table were in most cases found by this greedy algo- 
rithm. For n = 8, d = 5, the greedy algorithm gave a PA of size 26. However, 

P(8,5) > P(7,4) > 28 

by Theorem [6] Similarly, 

P(10, 7) > P(9, 6) > P(8, 5) > 28. 

Some other of the lower bounds are also determined using Theorem |6l They 
are marked by *. The upper bound is the Hamming type bound in Theorem [TTJ1 
or it's modified bound in Theorem 1111 Since P(n, 1) = nl for all n, this is not 
included in the table. 



Table 2: Bounds on P(n, d). 





d 


= 2 


d 


= 3 


d 


= 4 


n 


= dA 


- 1 




3 




3 




3 


n 


= rfH 


-2 


6 


- 24 




9 


9- 


- 12 


n 


= dA 


-3 


29 


- 120 


20 


- 34 


28 


-43 


n 


= dH 


-4 


90 


- 720 


84 


- 148 


68- 


- 166 


n 


= rfH 


-5 


582 


- 5040 


401 


- 733 


283- 


-4077 




d 


= 5 


d 


= 6 


d 


= 7 


n 


= d- 


h 1 




3 




3 




3 


7? 


= d- 


^2 


9 - 


- 12 


9- 


- 18 


9 - 


- 18 


n 


= rf + 3 


28* 


-43 


28* 


- 60 


28* 


-60 


n 


= d- 


h4 


95 - 


- 166 


95* 


- 216 


95* 


- 216 


n 


= d- 


h5 


236 


- 714 


236* 


-850 


236* 


-850 



5 Conclusion 

We give a number of constructions of permutations arrays under the Chebyshev 
distance, some with efficient decoding algorithms. We also consider an explicit 
mapping of vectors to permutations with efficient encoding/decoding. Finally, 
we give some bounds on the size of PAs under the Chebyshev distance. 
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